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ABSTRACT 

The Andromeda Galaxy (M31) is the closest large galaxy to the Milky Way, thus it is an important 
laboratory for studying massive dark objects in galactic halos (MACHOs) by gravitational microlens- 
ing. Such studies strongly complement the studies of the Milky Way halo using the the Large and 
Small Magellanic Clouds. We consider the possibilities for microlensing surveys of M31 using the next 
generation of wide field imaging telescopes with fields of view in the square degree range. We consider 
proposals for such imagers both on the ground and in space. For concreteness, we specialize to the 
SNAP proposal for a space telescope and the LSST proposal for a ground based telescope. We find 
that a modest space-based survey of 50 visits of one hour each is considerably better than current 
ground based surveys covering 5 years. Crucially, systematic effects can be considerably better con- 
trolled with a space telescope because of both the infrared sensitivity and the angular resolution. To 
be competitive, 8 meter class wide-field ground based imagers must take exposures of several hundred 
seconds with several day cadence. 

Subject headings: gravitational lensing — galaxies: halos — galaxies: individual (M31) — dark matter 



1. INTRODUCTION 

Almost twenty years ago Paczyhski (1986) proposed 
that dark objects in the Milky Way halo might be found 
by searching for the rare gravitational microlensing of 
stars in nearby galaxies: the Large and Small Magellanic 
Clouds (LMC and SMC) and also M31. This technique 
could thus be used to explore the dark matter problem in 
searching for these objects, dubbed MACHOs. Both pri- 
mordial black holes (a candidate for non-baryonic dark 
matter, known to dominate the clustering matter in the 
universe) and stellar remnants (a candidate for bary- 
onic dark matter, an unknown component accounting for 
roughly half of the baryons today) are interesting MA- 
CHO candidates. For a dark halo consisting of solar- 
mass MACHOs, of order one in one million LMC stars 
is being lensed (with a magnification of 30% or more) at 
any given time. These events are transient, developing 
over typically 100 days. A monitoring campaign could 
hope to detect these events, thus shedding light on the 
nature and composition of the Milky Way halo. Several 
groups have performed such campaigns over more than 
ten years, and have extended the technique to more dis- 
tant galaxies as well. 

The MACHO collaboration (Alcock et al. 2000) mon- 
itored the LMC for more than 8 years, concluding from 
the 5.7 year dataset that there is a population of lenses 
making a contribution of approximately 20% in objects 
of roughly half-solar mass. The EROS collaboration 
(Afonso et al. 2003) has monitored the LMC and SMC 
over a similar time period, and finds only an upper limit 
of 25% on the microlensing component. The MACHO 
collaboration result is intriguing in that while the possi- 
bility that the dark matter consists solely of these lensing 
objects is ruled out, there may be a subdominant non- 
baryonic dark matter component. Otherwise, the result 
speaks to the dark baryon problem. 

More than a decade ago, Crotts (1992) pointed out 
that the Andromeda Galaxy (M31) would be an excellent 
target for a microlensing survey, in particular because of 



the asymmetry expected in the rate of halo lens events 
due to the steep inclination of the M31 disk. Both the 
Milky Way and M31 halos could be studied in detail. 
Subsequently, this result was discussed in depth (Jctzcr 
1994; Baltz & Silk 2000; Kerins et al. 2001; Baltz, Gyuk 
& Crotts 2003). 

Microlensing surveys of M31 require image subtraction, 
as very few stars can be resolved from the ground. This 
is the regime of "pixel" lensing (Crotts 1992; Baillon et 
al. 1993; Tomaney & Crotts 1996; Gould 1996). Several 
collaborations, including MEGA (de Jong et al. 2004) 
(preceded by the VATT/Columbia survey (Crotts & 
Tomaney 1996, Uglesich 2001, Uglesich et al. 2004)), 
AGAPE (Ansari et al. 1999, Auriere et al. 2001, Calchi 
Novati et al. 2002, Paulin-Hcnriksson et al. 2003, Calchi 
Novati et al. 2003), and WeCAPP (Riffcser et al. 2003), 
have produced a number of microlensing event can- 
didates involving stars in M31. The results of the 
VATT/Columbia survey (Uglesich et al. 2004) indicate 
the presence of a lensing halo at the 20% level, excluding 
the case of only self lensing at greater than 95% confi- 
dence. 

Gould (1995) noted that the Hubble Space Telescope 
(HST) could perform a microlensing survey of the Virgo 
cluster using M87. Such a survey has been carried out 
(Baltz et al. 2004), placing only an upper limit that the 
Virgo halo is less than 60% lenses in the solar mass range. 

The current observational situation is uncertain. 
Ground based programs are moving forward, such as 
SuperMACHO, MEGA, and AGAPE. To understand 
populations of dark objects in galactic halos, wide field 
imaging telescopes currently being proposed and studied 
would have large advantages. In this paper we discuss the 
potential of wide field imaging for microlensing surveys 
of M31. For significant progress to be made, ground- 
based surveys would need to be performed with 8 me- 
ter class telescopes, or with night-long integrations on 
smaller telescopes. We discuss using the proposed LSST 
for this purpose. A space-based wide field survey offers 
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many advantages relative to the ground; we discuss us- 
ing the proposed SNAP space telescope for microlensing. 
Graff & Kim (2001) have previously discussed using the 
SNAP telescope for a microlensing survey of the LMC. 

2. MICROLENSING RATE 

Point gravitational lenses always exhibit two images of 
any source. The separation of these images for solar mass 
lenses within the Local Group tends to be microarcsec- 
onds, thus the term microlensing. Even with unobserv- 
able image splitting, the magnification can be large. For 
lenses and sources in relative motion, the magnification 
is transient, resulting in a microlensing event. In this sec- 
tion we describe the calculation of the rate of detectable 
microlensing events. 

2.1. Basics 

The lens equation can be written in terms of the Ein- 
stein radius and angle, which for a point mass lens are 
given by 
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where M is the lens mass, A is the distance to the lens, 
A is the distance to the source, and As is the distance 
between the lens and source. The lens equation relating 
source position 9 s and image position 9i is then given by 
(in angular coordinates with the lens at the origin) 



9 s = e I [i 
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This equation can be written with all angles written in 
units of the Einstein angle, in particular u — 9s/9e- The 
magnification was first given by Einstein (1936): 



A(u) = l + f(u% f(x) = 
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For ti«l, ArsI/b, thus large magnifications are pos- 
sible for point sources. 

A microlensing event has a characteristic timescale 
called the Einstein time, £e = Re/vj_, given in terms 
of the perpendicular velocity of the lens relative to the 
source vj_. For rectilinear motion, we find 
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where [3 is the minimum impact parameter of the lens, 
in units of the Einstein angle. For a star with unlensed 
flux A, the microlensing lightcurve is then 

F(t) = A + AiV(i)) = A + FJ(u(t)). (5) 

The Einstein time can be easily extracted from the 
lightcurve if A is known. The lightcurve shape contains 
some information on f3 (and thus A), though the de- 
pendence is weak, requiring very high quality data (e.g. 
Paulin-Henriksson et al. 2003). 

In the previous discussion we have assumed point 
sources. If the angular size of the source (9* = A/ A) is 
large, namely 9* > (39e, finite source size effects can be 
significant, as summarized by Yoo et al. (2004). 



2.2. "Pixel" Microlensing 

The unlensed flux of the source star is difficult to mea- 
sure. Significant blending is possible even for nearby 
sources (e.g. in the LMC). Most sources in M31 are highly 
blended, even with the resolution available from a space 
telescope. The full width at half maximum timescale, 
rather than the Einstein timescale, is the one usually 
measured. It is given by ty 2 = 2tEw((3), where 
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w(0) = WtfW 2 )) - P 2 , 

with limiting behavior 

w(j3 < 1) = 0y/3, w(j3 > 1) = (3\f V2-1. (7) 

We find that ty 2 ~ £e P for ah values of f3. 

The degenerate form (high magnification limit /? — ► 0) 
of the microlensing lightcurve is 
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where B is the baseline flux and AAax = A<5(/3) ~ 
A//3 is a fit parameter expressing the maximum increase 
in flux from the lensed star. In the absence of blending 
B = A- The "pixel" lensing regime in where sources 
are only resolved when they are lensed. The measured 
parameters are AAiax and ty 2 . Without knowing A, 
measuring the Einstein time £e requires measuring /?, 
which provides only a small correction to the degener- 
ate lightcurve. High quality data is required (Paulin- 
Henriksson et al. 2003). We proceed using the distribu- 
tion of microlensing events in ty 2 , though not as useful 
as that in £e- We note that AAax^i/2 ~ A^e in the 
high magnification regime, so Einstein times can be re- 
covered statistically, given that the distribution of A is 
known (Gondolo 1999). This distribution in A is sim- 
ply the stellar luminosity function, and at least the first 
non-trivial moment is known: the surface brightness fluc- 
tuation (SBF) magnitude. 

In practice the lightcurves are extracted by image dif- 
ferencing. A reference image is constructed from some 
number of individual exposures, and is then subtracted 
from the individual images. In this way the baseline is 
removed, and the only remaining sources are variable. 

2.3. Microlensing Rate 

The timescale distribution for microlensing events is 
expressed as the integral along the the line of sight of the 
lens density times the cross section (Griest 1991; Baltz 
& Silk 2000): 
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where A = xD s , uj-e = A/(w c £e), ^0 is a modified 
Bessel function of the first kind, v c is the circular ve- 
locity of the lens population (v c — 0"v/2), and r\ — Vt/v c 
is the transverse velocity of the line of sight relative to 
the lens population. Lenses are assumed to be drawn 
from a maxwellian velocity distribution. The transverse 
velocity due to the motion of source and observer, as- 
suming that the source velocity is prescribed and not 
taken from a distribution (e.g. disk sources whose ran- 
dom velocities are much less than their bulk velocities 
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in the halo frame). This equation does not account for 
detection efficiency, which can be expressed as a prob- 
ability to detect an event with certain fit parameters 
P[ti/2,Qmax,0\, where Q max oc A.F max is just the maxi- 
mum flux expressed as a signal-to-noise ratio. The rate 
of detectable events is then just 

1T= dp^P[2tvw(p),QJ((3),0\. (10) 

Note that the fit parameters ty 2 and AF max have been 
replaced with more physical ones depending on [3. Here, 
Q-,, (oc F+) is the naive (photon counting) significance 
with which the source star is detected if blending is ig- 
nored. Since dTo/dt E is independent of f3, the rate can 
be written as 

-^=/3eff(*E,Q*) ^(M lens ). (11) 

The distribution is integrated over £g, and then over the 
mass function of lenses and the luminosity function of 
sources to arrive at the total observable rate. Note that 
t~E may be unknown event by event, but integrating over 
all £e yields the total observed event rate. The distribu- 
tion dT/dti/2 can be derived similarly, 

dT f°° dV f fya \ P[t V 2,QJ{P),0\ 
dt l/2 J P dt E \2w(0)J 2w((3) 

(12) 

We include finite source effects using a simple prescrip- 
tion (Baltz & Silk 2000). For a given source, we can 
determine the required magnification to give a lensed 
flux large enough to pass a detection criterion. For a 
given source diameter the maximum magnification 
is a function of x; as x — > 1, A max — > 0. We can solve 
for the largest x allowing the required magnification, and 
truncate the x integral in equationElaccordingly. A quick 
estimate shows that at most a handful of events might 
be observed toward M31 with significant finite source ef- 
fects. 

We have assumed that the source velocities have 
no random component. If the source velocities are 
maxwellian (as should be approximately the case for a 
bulge), the velocity integrals can be separated out. The 
outcome is the same if the identification w^(lens) — > 
u^(lens) + x 2 v1 (source) is made. 

2.4. Monte Carlo Event Rates 

For a given microlensing survey, we can perform a 
Monte Carlo simulation to determine the event detec- 
tion probability P\ty 21 Qmaxi P]- For this work it will be 
safe to assume that the Poisson errors in photon count- 
ing will be approximately gaussian, i.e. the number of 
photons expected is significantly larger than 1. 

The time sampling in a given survey is used to gen- 
erate simulated microlensing lightcurves with fit param- 
eters ty2, Qmax, P- Random peak times to are chosen. 
Data points with Poisson errors are generated for each 
time sample. The generated lightcurves are then passed 
through a set of criteria for detection. In this way, the 
probability P is calculated for a set of fit parameters. 
This process is repeated, generating a table of probabil- 
ities as a function of the fit parameters, to be used in 
eauationllOl 



In this work we use the following detection threshold. 
A total of lOcr is required, collected over N consecutive 
samples that are 10a /VlV above the baseline. We let 
N = 2, 3, 4, 5, 8, allowing for short bright events as well 
as longer dimmer events. 

3. MODELING OF M31 

We use the microlensing model of M31 described in 
Baltz, Gyuk & Crotts (2003). This model consists of 
sources in the bulge and disk of M31, and lenses in the 
bulge, disk, and halo of M31, as well as Milky Way halo 
lenses. 

3.1. Galactic Halos 

We model the dark halos of the Milky Way and M31 as 
cored isothermal spheres, having asymptotic 1/r 2 density 
profiles, giving flat rotation curves. Lenses making up a 
fraction /mw,M31 of the dark halo have a space density 
given by 

v 2 yi - ^/ (q sin- 1 yi - ^) 

p{x, y,z) — f — £ — — — , (13) 

iirG x z + y z + (z/qy + r£ 

where v c is the asymptotic rotation speed, r c is the core 
radius, and q is the flattening, with q = 1 indicating no 
flattening. We assume that v c = 220 km s -1 for the 
Milky Way halo and v c = 240 km s" 1 for the M31 halo. 
We vary the core radius for M31 as r c =0.1, 1, 2, 5, 10 
kpc, and take flattening values q — 1 and q = 0.3. We 
assume that the halo velocity distribution is maxwellian. 
We impose halo cutoffs at a distance of 200 kpc from the 
respective centers. 

3.2. Stellar Populations 

Sources in M31 are taken to reside in a two component 
model consisting of a double-exponential disk (Walter- 
bos & Kennicutt 1988; Gould 1994) and a bulge (Kent 
1989). Details of this model have been discussed pre- 
viously (Baltz, Gyuk & Crotts 2003). Self lensing be- 
tween these components is computed (all cases: bulge- 
bulge, bulge-disk, disk-bulge, and disk-disk), serving as 
the background for the detection of a lensing halo. The 
M31 disk and bulge are then considered as sources for 
halo lensing events. 

3.3. Mass and Luminosity Functions 

We adopt the Chabrier (2001) mass function for stellar 
lenses in M31, both in the disk and the bulge. For halo 
lenses we assume a delta function mass function, though 
varied between 10 _3 M Q and 10 Mq. For the luminosity 
function of stars, we consider the bulge and disk sepa- 
rately. For the disk, we use the i?-band (for ground based 
surveys) and J-band (for space based surveys) luminos- 
ity functions of Mamon & Soneira (1982). The bulge is 
more complicated, as data are scarce. We take an average 
(in log dN/dM) of F-band and /-band data for M R > 
(using Terndrup, Frogel & Whitford (1990) for /-band 
and Holtzman et al. (1998) for y-band, both toward 
the Milky Way bulge). For Mr < we attach a power 
law slope of 0.59 taken from the MACHO project data 
(Alves 2001). For the J-band, we use an isochrone for an 
old metal-poor population (Lauer & Mighell 2003) with 
surface brightness fluctuation magnitude Mj — —3.43, 
which is adequate for this study. 
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The full dataset of a microlcnsing survey would allow 
for better measurements of the stellar luminosity func- 
tions, at least in terms of the SBF magnitude. A space 
based survey would obviously be considerably better for 
this purpose than a ground based one, simply due to the 
considerably better resolution. 

4. SAMPLE WIDE FIELD SURVEYS 

From the microlensing model of M31, we now proceed 
to study wide field surveys. Some experimental parame- 
ters arc needed, namely the sensitivity and angular reso- 
lution of the imaging systems in question. In the follow- 
ing, we discuss a ground based plan and a space based 
plan. 

4.1. Ground Based Survey with LSST 

The Large Synoptic Survey Telescope (LSST) is a pro- 
posed instrument for mapping the entire accessible sky 
every several days (Tyson et al. 2004). As such, it will 
discover microlensing in M31 and elsewhere. We assume 
some basic experimental parameters as follows. The 
LSST focal plane is very large, of order 10 square de- 
grees. All of the interesting region of M31 is easily cov- 
ered in one field. One caveat is that the site of LSST 
is as yet undecided, and the visibility of M31 from the 
southern hemisphere is limited. However, the discussion 
carries over to any wide-field imaging system with expo- 
sure times scaled appropriately. 

The LSST aperture is an annulus, with an effective di- 
ameter of 6.5 meters. One possible plan for the LSST 
survey is to scan the entire visible sky in 4 colors every 4 
days, with 10 second exposures, possibly in pairs. We will 
assume an i?-band photometric zero point of rnjf = 27. 4 
mag (the flux giving one photoelectron per second), im- 
plying a 10 second exposure zero point of m° R = 29.9 
mag, giving one photoelectron over the exposure. Fur- 
thermore, we assume that the nose level is twice the pho- 
ton counting noise, as is typically achievable in ground 
based microlensing of M31 (the degradation is possibly 
due to unresolved variable stars). We use a cadence of 
4 days to do the Monte Carlo simulation described pre- 
viously. This is an oversimplification, as bad weather is 
not taken into account, but acceptable for our purposes. 

The resolution can be characterized by one number, 
usually taken to be the seeing. More useful for mi- 
crolensing is SIpsf = where tpi is the normal- 
ized PSF (Gould 1996). For a gaussian PSF with seeing 
(full width at half maximum) do, it is easy to show that 
^psf = In 4. We assume that the typical seeing at 
the LSST site (still being discussed) will be 0.75", giving 
fipgp = 1.27 square arcseconds. 

We take the surface brightness /jr from the M31 model 
(Baltz, Gyuk & Crotts 2003). The background sky 
brightness is assumed to be uniform, with /i s k y = 21 
mag arcsec -2 in the i?-band. Lastly, we assume that 
the distance modulus to M31 is D = 24.5, corresponding 
to a distance of 795 kpc. We can now express the con- 
version to signal-to-noise ratio as follows, for a star with 
apparent magnitude wlr magnified by a factor A = 1 + 5, 

1 iQ-0A(m a -m° a /2) g 

Q ~2 v /n P 7 F (lO- - 4 "* + io-o-4^ ky ) ' ^ 

and the factor of 1/2 is due to the fact that the noise 
level is twice the photon counting noise, usually the case 



for observing M31 from the ground. This is likely due to 
unresolved variable stars, e.g. RR Lyrae (Crotts 2004). 
This expression assumes that the noise in the reference 
image used in image subtraction is much less than the 
noise in individual images, as is usually the case. 

4.2. Space Based Survey with SNAP 

The SuperNova Acceleration Probe (SNAP) is a pro- 
posed wide field space telescope chiefly aimed at mea- 
suring the acceleration of the universe by Type la super- 
novae and by weak gravitational lensing (Aldcring et al. 
2004). 

The SNAP aperture is planned to be 2 meters in diam- 
eter. The focal plane covers 0.7 square degrees, though in 
a strange shape. There are 9 colors, with filters mounted 
directly on the focal plane. There are 4 optical and 4 in- 
frared sections. The optical sections are broken up into 
a 6x6 grid, with 6 filters, while the infrared sections are 
broken up into a 3x3 grid, with 3 filters. The current plan 
is for the 9 filters to start at B-band (central wavelength 
4400 A), rcdshifting by a factor of 1.15 going redward, 
the longest filter having a central wavelength of 1.346 
fim. This range covers the usual color system BVRIJ 
and slightly redder. Microlcnsing rates tend to be largest 
in the reddest filters, thus we will assume that J-band 
is used for the primary event detection. Converting the 
quoted sensitivity (Aldering et al. 2004) from AB magni- 
tudes, we use 77ij P = 23.7 mag, and for a 100 second ex- 
posure (the lower limit based on bandwidth) m°j = 28.7 
mag. We conservatively estimate that fipsF = 0.1 square 
arcseconds. 

A microlensing survey with SNAP would take time 
away from the main surveys, but the resources required 
are modest. To cover the central region of M31, 27 point- 
ings would be adequate. This would cover each position 
in the 3 infrared filters and 3 of the optical filters. The 
optical coverage could be staggered, so that in two visits, 
all points would be covered in the 3 infrared filters twice 
and the 6 optical filters once. Bandwidth to the tele- 
scope limits the exposure times to be more than about 
100 seconds, so we assume 100 second exposures. This 
means a visit would consist of 2700 seconds of integra- 
tion, plus overheads which may be as large as 50%. M31 
is at low ecliptic latitude (+33.3°), thus even from space 
it is only visible for part of the year. 1 Constraints from 
the roll angle of the telescope would likely limit the vis- 
ibility to 2 periods of 90 days in a year, when the sun 
is 90° away from M31. We consider surveys of 90 days, 
with cadences of 2 days and 4 days (in the optical filters 
4 days and 8 days). 

We adapt the M31 model to J-band simply estimat- 
ing that R — J = 1.67 for the disk and R - J = 1.72 
for the bulge. We assume a sky brightness fij = 21.7 
mag arcsec -2 . Because many stars in M31 are resolved 
at this resolution, we must be careful in estimating the 
noise level. We take a conservative approach, adding the 
star brightness to the background. A more correct ap- 
proach would divide the surface brightness into resolved 
and unresolved components, thus there would be no dou- 
ble counting. We connect to signal-to-noise ratio simi- 

1 As pointed out by Graff & Kim (2001), the LMC is an at- 
tractive target for a microlensing survey with SNAP for just this 
reason - it lies within 5° of the south ecliptic pole. 
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TABLE 1 
Expected Number of Events 



survey 


self 


MW 


M31 


self + 20% halos 


SNAP 2 day 


165 


240 


635 


340 


SNAP 4 day 


110 


185 


470 


240 


LSST 10 s 


18 


33 


79 


10 


LSST 100 s 


99 


155 


400 


210 


LSST 1000 s 


465 


590 


1590 


900 



Note. — Expected number of events in mock surveys. We 
assume a core radius for the halos r c = 2 kpc and no flattening 
(q = 1). The MACHO mass is fixed at -1/2 dex solar (.32 M©). 
Numbers for the individual halos assume a 100% MACHO fraction. 
The SNAP survey covers 90 days, while the LSST numbers are 
given per full year of observations - for M31 this would likely take 
3 years. In addition to the baseline 10 second exposure for an 
LSST survey, we also give the rates for identical surveys (4 day 
cadence) taking 100 second and 1000 second exposures as well. 
The LSST numbers can easily be scaled to different telescopes, 
e.g. an unobscured 3.2 meter telescope has 1/4 the area, thus 
the numbers are appropriate for 40 second, 400 second, and 4000 
second exposures. 



larly to the ground-based case, 

1Q -0A(mj-m O j/2) r 

Q= — (15) 

Am j + Q psF (10-0.4 MJ + 10 -0.4 Msky ) 

The resolution of a space-based telescope such as 
SNAP provides an important advantage to a microlens- 
ing survey of M31 in that some of the stars will be re- 
solved. For such stars, the Einstein timescale can be 
measured directly. We use a conservative criterion for 
determining which stars are resolved: if a star is 10 times 
brighter than the surface brightness contained in f2psF, 
then it is resolved. This can be expressed as 

m < — 2.5 — 2.5 log f2psF — ► m < Mi (16) 

where m is the apparent magnitude of the source, and 
fi is the surface brightness of the galaxy (in the same 
band). 

4.3. Number of Events in Mock Surveys 

For the sample surveys in both LSST and SNAP we 
provide the expected number of events. These are listed 
in Tabled We compute rates for flattenings q — 1 and 
q = 0.3, and for halo core radii r c =0.1, 1, 2, 5, 10 kpc. 
We use a delta function for the halo lens mass function, 
varying between 1O _3 M0 and IOMq, in intervals of 1/8 
dex. Table n illustrates the fiducial case q = 1, r c = 2 
kpc, M = —1/2 dex solar. We also compute the distri- 
bution in timescales for these events, both dT/dt-E and 
dT /dt i/2- In the case of SNAP, we separately calculate 
the number of events with £e determination according to 
equation ^] and those with only ty 2 determined. 

5. MAXIMUM LIKELIHOOD 

Having calculated the expected number of microlensing 
events in various surveys, we would like to understand 
how well the parameters of the microlensing model can 
be measured. We use a maximum likelihood technique 
first developed for the VATT/Columbia survey (Uglesich 
et al. 2004), as we describe below. 



5.1. Simulated Datasets 

We simulate surveys of M31 based on the expected 
event rate as calculated in the previous section. The 
event rates are binned in position and log timescale with 
intervals of 1 arcminute and 1/8 decade respectively, for 
each component of the total lensing rate. For a given set 
of parameters, such as lens mass and halo lens fraction, 
a total expected lensing rate Tj is constructed in each bin 
across the field of view and in timescale. Simple Poisson 
statistics are then used to generate a number of events in 
each bin (usually zero as there are many more bins than 
expected events). The simulated data is then a set of 
number of events in bins, {n^}. It is simple to generate 
many realizations for each set of parameters. In fact 
we generate two parallel sets of bins, for the expected 
event rate with and without a measurement of £e- In 
Figure^we illustrate one such realization, for a halo lens 
fraction of 20% and lens mass —1/2 dex solar (0.316M©) 
. We assume the 2 day cadence SNAP survey of 90 day 
duration, and show the size of SNAP focal plane in the 
infrared filters (the optical filters are the mirror image). 
Note that we have separated source stars into two groups 
according to eauation ll6l and then separately generated 
realizations for events with measured £e and those with 
only ty2 measured. In Figure [21 we show this simulated 
dataset in more detail, both illustrating the timescales 
of events, and also focusing on the central region, which 
has the highest event density, and where star-star lensing 
dominates. 

5.2. Likelihood Function 

The likelihood function is clearly just the product of a 
Poisson distribution for each bin, for a number of events 
ni, where the expectation is Tj, 

ni 

i 

The expected number of events in each bin is a function 
of the parameters of the lensing model. The most impor- 
tant are lens mass and halo lens fraction. The likelihood 
function is maximized over these parameters. We can 
define \ 2 — —2 ln£, though we don't expect this quan- 
tity to be distributed in the usual way as this regime of 
Poisson statistics is far from gaussian, 

x 2 = 2 J2 in - n i lnT * + M^O] ■ (18) 

i 

For each set of input parameters (lens mass, lens frac- 
tion, core radius, flattening), we compute the likeli- 
hood function C. Marginalizing over halo shape (core 
radius and flattening), we construct contours of equal 
likelihood. We integrate the likelihood function over 
the lens mass and halo lens fraction to identify the 
contours containing probabilities corresponding to the 
usual la (68.269%), 2a (95.450%), 3ct (99.730%), and 
4a (99.994%) confidence levels. We also marginalize over 
each axis in turn to arrive at the one-dimensional confi- 
dence intervals for lens mass and halo lens fraction. 

In Figure [3] we illustrate the likelihood contours for 
the cases of 20%, 5%, 2% and 0% lensing halos. For each 
case, 10 realizations were generated, and likelihoods av- 
eraged with equal weight for each. In the case of 2%, 
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Fig. 1. — SNAP focal plane with simulated dataset, assuming a halo lens fraction of 20% with mass —1/2 dex solar. The cadence is 2 
days, with a total survey duration of 90 days. The positions of the infrared detectors are indicated, showing the scan across the field. The 
infrared detectors are individually 1/3 the width of the illustrated squares — each square thus requires 3 exposures in the scan. With the 
edge effects, 26 exposures are required to cover the central 4x5 chip region, which is boxed. To stagger the position in alternating scans for 
the optical filters (whose positions are the mirror image of those shown), one additional exposure is required, making 27. Dots mark the 
simulated events, and diamonds mark those few simulated events passing the 4e measurement criterion. Note the considerable asymmetry 
between the near (bottom) and far (top) sides. 



the likelihood contours vary significantly between real- 
izations, indicating that this is about the limit where a 
lensing halo could be detected. In Figure 0] we illustrate 
several of the separate realizations for the 2% case. 

Finally we consider the marginalized likelihoods for 
lens mass and halo lens fraction. We can simply integrate 
the 2-dimensional likelihood plots along each dimension 
in turn, and find the confidence intervals. These intervals 
are given in Table [3 

6. CONTROL OF SYSTEMATICS 

Wide field field microlensing surveys of M31 offer the 
possibility of controlling systematic difficulties consider- 
ably more robustly than has been possible to date. A 
space based wide field survey would be especially well 
suited to this task. 

The relevant systematics in a microlensing survey all 
involve backgrounds. The two most important are the 
problems of variable star rejection and self lensing. Some 
classes of variable stars, in particular the Mira variables, 
can exhibit lightcurves that mimic microlensing, and re- 
peat only after periods of months or years. Identifying 
these classes of variable stars is thus crucial in that false 
events must be avoided. Self lensing (lensing of stars by 
stars) is a significant background for microlensing sur- 
veys. In the mock surveys discussed here, in fact half of 
the events are self lensing assuming 20% lensing halos. 
The self lensing can be separated from the halo lensing 



by its spatial distribution, but this must be understood 
well. 

General purpose multicolor wide field imagers have sig- 
nificant power in understanding and accounting for these 
systematics. Wide field imaging makes for a highly ho- 
mogeneous dataset. In particular for M31, measuring 
the near / far asymmetry is complicated by non-uniform 
time sampling, and due to observational difficulties, the 
sampling of near and far side fields is rarely even the 
same. In practice this can be modeled, but this certainly 
introduces a systematic error. In addition, the presence 
of dust lanes spoils the near-far symmetry of stellar back- 
grounds. Multicolor observations of the entire field, es- 
pecially into the infrared, will allow this effect to be ac- 
counted for much more robustly. In fact, such observa- 
tions, especially with space-based resolution, will allow a 
more accurate self lensing model to be constructed. 

Having multiple colors, especially into the infrared, 
will also allow very sophisticated identification of vari- 
able stars. Most simply, the known classes of trouble- 
some variable stars tend to be red; having multiple col- 
ors for each candidate event, including infrared colors, 
will allow identification of likely variable stars with some 
confidence. In particular, Miras and semiregulars (both 
classes are red supergiants) are clearly identifiable in J- 
band (accessible to a space-based survey). Equally as 
important, following a candidate microlensing event in 
multiple colors allows the crucial test of gravitational 
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Fig. 2. — Detail of simulated dataset shown in Figure^ along with contours of expected event density. Again, the halo lens fraction is 
20% and the lens mass is —1/2 dex solar. The bottom panel shows the full simulated region, while the top panel shows the M31 bulge 
region. Circles indicate simulated events, with timescales in days of (smallest to largest, with colors online only) 0.3d < ty2 < Id ( rec 0, 
Id < ty 2 < 3d (orange), 3d < ty % < lOd (yellow), lOd < t 1/2 < 30d (green), 30d < t 1/2 < lOOd (blue), lOOd < ty 2 < 300d (purple). 
Diamonds indicate events where could be measured. The contours indicate expected events per square arcminute, and are spaced by 
1/2 decade (so "3" means vTO in the contour labels). Again, the near-far asymmetry is clear from the contours. 



lensing, that the color remains constant in time. Vari- 
able star outbursts usually have color variations, which 
would rule out the possibility of microlensing. 

A significant amount of baseline is required to identify 
long period variable stars. In a ground based survey like 
LSST, this will not be an issue, as at least short exposure 
data will be taken for the duration of the survey, likely 
several years. A space based survey would likely need at 



least some ground based support. While it would be pos- 
sible to take several exposures in each visibility period, 
more complete coverage would be desirable. Judging by 
event rates, infrequent exposures of several hundred sec- 
onds on an 8 meter class telescope would have adequate 
sensitivity to search for additional peaks in a candidate 
microlensing lightcurve. 
In the large surveys discussed in this paper, a signif- 




Fig. 3. — Likelihood contours for several halo fractions. Top left: 20%, Top right: 5%, Bottom left: 2%, Bottom right: 0%. The contours 
correspond to the usual la, 2a, 3a, 4<r confidence levels. In each case, 10 datasets were generated and the likelihoods were calculated. 
Then the results were averaged. The average results are quite representative in each case except for the 2% halo, where the contours for 
specific realization vary significantly — this is an indication that a 2% halo is about the level of sensitivity for this program. 



icant number of events are expected from stellar lenses 
alone. A considerable fraction (of order 10%) of these 
would be expected to display the caustic crossing signa- 
ture of a binary lens (Mao & Paczyhski 1991; Baltz & 
Gondolo 2001). Caustic crossing events exhibit unique 
lightcurves, thus they can be identified as microlensing 
with high confidence. The number of binary events can 
then be used to estimate the expected number of single 
events related to the same lens population. This can be 
used as a calibration of the self lensing model. 




Fig. 4. — Likelihood contours for a 2% lensing halo in several realizations. In four cases, there is a solid detection of a lensing halo, in 
two others, a weak detection, and the remaining two, no detection at all. A solid detection was thus found in half of the realizations. For 
this reason we argue that a 2% halo is roughly the limit that the SNAP survey could detect. 



TABLE 2 

Confidence Intervals 



lens fraction 


lens fraction 


lens fraction 


lens mass 


lens mass 


lens mass 


lens mass 


(percent) 


(percent) 


(percent) 


(solar masses) 


(solar masses) 


(dex solar) 


(dex solar) 


input 


la 


2a 


la 


2a 


la 


2a 


20 


16 — 23 


13 — 28 


0.19 — 0.35 


0.14 — 0.49 


-0.72 0.46 


-0.84 0.31 


5 


3 — 6 


2 — 8 


0.19 — 0.45 


0.13 — 0.73 


-0.73 0.34 


-0.89 0.14 


2 


0.8 — 3.8 


0.0 — 5.8 


0.14 — 0.51 


0.079 — 2.6 


-0.85 0.29 


-1.10 1-0.42 





0.0 — 1.0 


0.0 — 2.2 


0.017 — 2.1 


0.0013 — 7.8 


-1.76 h0.31 


-2.88 hO.89 



Note. — Confidence intervals for measuring the lens mass and halo fraction in the SNAP survey with 2 day cadence discussed in the 
text. On average, the 2% lensing halo is not detected with high confidence, though in several cases it was (see FigureEJ. The 20% halo was 
recovered with errors of 3.5% in halo fraction and 0.13 dex in lens mass. The 5% halo was recovered with errors of 1.5% in halo fraction 
and 0.20 dex in lens mass. In the case of no lensing halo, the 95% limit on the halo lens fraction was put at 2.2%, a very considerable 
improvement on the capabilities available now or in the near future. 
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7. DISCUSSION 

We have calculated the sensitivity to microlensing of 
two proposed wide field imaging systems, and found that 
even with modest programs these would eclipse anything 
currently available. More importantly, such wide field 
multi color surveys would address the systematic diffi- 
culties with convincingly identifying microlensing more 
robustly than has been possible to date. A space based 
survey such as SNAP would be especially powerful for 
this purpose due to the availability of high resolution 
and infrared imaging. 

The chief science goal of the microlensing surveys de- 
scribed in this paper would be to understand the popula- 
tions of stellar-mass objects (if any) associated with the 
halos of spiral galaxies. This speaks to the issue of the 
dark baryons, still unidentified, and making up roughly 
half of the baryon budget in the present-day universe. 

In addition to dark halo objects, other science products 
of such surveys would include a census of variable stars 
(in several filters), and other transient events. One in- 
teresting possibility is the detection of binary lens events 
with extreme mass ratios, indicative of giant planets. 
This has been recently achieved toward the galactic bulge 
with a solid detection by the OGLE and MOA teams 
(Bond et al. 2004), followed by another detection by 



OGLE (Jaroszynski et al. 2004). Detecting planets would 
necessarily involve reducing the survey data in real time, 
so that microlensing alerts could be issued, and candi- 
dates could be followed closely, either from the ground 
or in space. 

We conclude this paper with a simple observation. The 
wide field (of order square degree) imaging systems being 
discussed now have obvious applications for cosmological 
surveys. In contrast, there are only a handful of objects 
on the sky with angular sizes commensurate with these 
imagers, including the Milky Way bulge and disk, the 
Magellanic Clouds, and M31. It seems warranted that at 
the very least, modest studies of these very few nearby 
objects should be undertaken with the upcoming wide 
field technologies. 
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